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Abstract. We discuss an alternative approach to quintessence modifying the usual 
equation of state of the cosmological fluid in order to see if going further than the 
approximation of perfect fluid allows to better reproduce the available data. We 
consider a cosmological model comprising only two fluids, namely baryons (modelled 
as dust) and dark matter with a Van der Waals equation of state. First, the general 
features of the model are presented and then the evolution of the energy density, the 
Hubble parameter and the scale factor are determined showing that it is possible to 
obtain accelerated expansion choosing suitably the model parameters. We use the the 
data on the dimensionless coordinate distances to Type la supernovae and distant radio 
galaxies to see whether Van der Waals quintessence is viable to explain dark energy 
and to constrain its parameters. We then compare the model predictions with the 
estimated age of the universe and the position of the first three peaks of the anisotropy 
spectrum of the cosmic microwave background radiation. 



PACS numbers: 98.80.-k, 98.80.Es, 97.60.Bw, 98.70.Dk 
1. Introduction 

In the last few years an increasing bulk of data have been accumulated favouring the 
scenario of a spatially flat universe dominated by some form of dark energy. A first 
strong evidence came from the Hubble diagram of type la supernovae (hereafter SNela) 
that turned out to be best fitted by spatially flat accelerating cosmological models 
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including a non trivial component pQ . On the other hand, the results from the observed 
first and second peak in the cosmic microwave background radiation (hereafter CMBR) 
spectrum strongly suggested that the geometry of the universe is spatially flat [21 IS]- 
When combined with the data on the matter density parameter VLmi these results lead 
to the conclusion that the contribution Qx of the dark energy is the dominant one, 
being (f2 M , fl x ) — (0.3, 0.7). This picture of the universe has been further strengthened 
by an increasing sample of high redshift SNela [3J |3] and most precise and extended 
measurements of the CMBR spectrum [HI Ej . 

After the discovery of these evidences of a spatially flat and accelerating universe, an 
overwhelming flood of papers, presenting a great variety of models for the explanation of 
this phenomenon, has appeared. The simplest explanation is the cosmological constant 
jHj which is able to fit the SNela data with good confidence, but it is also plagued by 
many problems on different scales. This situation has strongly encouraged the search for 
alternative approaches which now ranges from minimal coupled scalar fields, to strings 
and anthropic principle (see, e.g., |0] and references therein). 

In a recent work (TJIj, some of us have introduced a new approach to the problem 
modifying the usual equation of state of the cosmological fluid in order to see whether 
going further than the standard approximation of cosmological perfect fluid allows to 
reproduce the available data. The authors have considered a standard cosmology with a 
Van der Waals equation of state for matter without any other kind of energy source. This 
approach is extremely "natural" and "obvious" since it starts from the consideration 
that the universe, in its evolution, is not always well described by a perfect fluid in 
the forms of radiation or non- interacting dust. In our opinion, before adding an exotic 
and mysterious dark energy into the cosmic pie, it is worth wondering whether the 
observed acceleration of the universe can be implemented using the minimal number of 
fluids§. The price to pay is taking into account a more complicated equation of state and 
motivate it physically. From elementary thermodynamics, we know that a real fluid is 
never perfect [TI]. Moreover, it is also well known that the perfect fluid equation of state 
p = 7p with 7 = is just a rough approximation of cosmic epochs capable of describing 
stationary situations where phase transitions (e.g., from radiation dominated to dust 
dominated regions) are not considered ^2j. On the other hand, the only thing we know 
about dark energy is that it gives rise to an accelerated expansion, but there are no 
hints about its nature. Due to this situation, it is worth asking whether dark energy is 
indeed needed or, on the contrary, the observed acceleration is driven by standard dark 
matter provided that its equation of state is more realistically treated. 

A Van der Waals fluid could be a first step toward the goal to get a whole dynamics 
where 1.) only observed fluids are taken into account, 2.) phase transitions occur in the 
framework of the same evolution, 3.) accelerated and decelerated periods depend on 
the relative values of the parameters of the state equation with respect to the pressure 
and matter energy density which are functions of time. It is worth noting that a similar 

§ Following Newton "Hypotheses non jingo ". 
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approach has yet been explored in ^3] where the author considers a mixture of two 
fluids, using the perfect gas equation of state for the matter and the Van der Waals one 
for the dark energy. However, our approach is radically different since the model we 
consider is made out of matter only. Is is worth noting that the term matter usually 
refers to both baryons and dark matter and both these substances are described by 
the same equation of state. Actually, we have a direct knowledge of the properties 
of baryons only and indeed their equation of state is well described (on cosmological 
scales) by the dust approximation. On the opposite, the nature of the dark matter is still 
completely unknown so that, a priori, nothing prevents us from exploring the possibility 
that its properties call for a more general equation of state such as the Van der Waals 
one. The aim of the present paper is to explore further this approach in order to see 
whether a cosmological fluid with Van der Waals state equation can be reconciled with 
observations. We have thus constrained the effective parameters of the theory looking 
for cosmological models which: i) can admit a nowaday accelerating universe; ii) satisfy 
the constraints on the estimated age of the Universe, Hi) are able to fit the data on the 
dimensionless coordinate distance to SNela and radio galaxies. 

The paper is organized as follows. In Sect. 2, we briefly summarize the main features 
of Van der Waals equation of state and choose the set of parameters which are best suited 
to assign the model. Dynamics of the model is described in Sect. 3 where we determine 
the evolution of the energy density, the Hubble parameter and the scale factor. Matching 
with observations is performed in Sect. 4 where the data on the dimensionless coordinate 
distance to SNela and radio galaxies are used to select among the models. In Sect. 5, we 
compare the model prediction for the age of the universe with the estimated one, while, 
in Sect. 6, we evaluate the position of the first three peaks of the CMBR anisotropy 
spectrum and qualitatively speculate on how structure formation could take place in 
the Van der Waals quintessence scenario. Sect. 7 is then devoted to the discussion of 
the results and conclusions. 

2. The Van der Waals equation of state 

The dynamical system describing a Friedmann-Robert son- Walker (FRW) cosmology is 
given by the Friedmann equations : 



a 



AixG 



(Pb + Pdm + 3pdm) , 



(1) 



a 



3 




(2) 



and the continuity equations for each of the two fluids : 



Pi + 3H (pi+Pi) = , 



(3) 



where H = a/a is the Hubble parameter, the dot denotes the derivative with respect to 
cosmic time and k = —1, 0, 1 is the spatial curvature constant respectively for open, flat 
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and closed universes. Eqs.((T]), (J2J) and are derived by the Einstein field equations and 
the contracted Bianchi identities || assuming that the source of the gravitational field is 
a a mixture of baryons with energy density pb and pressure pb = and dark matter with 
energy density pbm and pressure Pdm- To close the system and determine the evolution 
of the scale factor a and of the other quantities of interest, the equation of state of the 
dark matter fluid (i.e. a relation between p DM and Pdm) is needed. 

In the standard cosmology, one assumes that the dark matter may be described 
as a perfect fluid so that the equation of state is pdm = IPdm where < 7 < 1 is 
the so called Zel'dovich interval. However, this is only an approximation which is not 
always valid and which does not describe the phase transitions between the successive 
thermodynamic state of the cosmic fluid. Several times, for instance at equivalence, 
two phases had to exist together. In these simple description by a perfect 

fluid equation of state is not realistic. An immediate generalization can be achieved by 
taking into account the Van der Waals equation of state which describes a two phase 
fluid. Also in this case, we have an approximation, but the consequences on dynamics 
are interesting [TU1 HB] . Hence we assume that the equation of state is : 

PvdW = 1 _ p pvdW ~ a Pvdw > ( 4 ) 

which reduces to the perfect fluid case in the limit a, (5 — > 0. Hereafter, we will denote 
with the subscript b (VdW) all the quantities referring to the baryons (the Van der 
Waals dark matter). In standard units, the a and (3 coefficients may be rewritten as: 

a = ZpcPc 2 ; P = (3p c ) _1 , (5) 

where p c and p c are the density and the pressure of the cosmic fluid at the Van der 
Waals critical point. The critical values are the indications that the cosmic fluid 
changes its phase at certain thermodynamic conditions. As a consequence, the three 
parameters (p c ,Pc,j) are n °t independent from each other. Inserting Eq.(j3J) into Eq.Q 
and considering the situation at the critical point, one gets: 

Pc = ^IPc , (6) 

so that the number of independent parameters is now reduced to two, which are (p c , 7). 
Using Eq. © we may rewrite Eq. (j3J) as : 

P 3777 9 2 

— = o 7 

Pc 3-77 8 

having defined the dimensionless energy density 77 = pvdw/Pc- It is worth noting that : 

9 

Pvdw ~ IPvdw for r? << 1 , p V dw ~ -^IPvdw for V » 1 • 

|| Actually, the Bianchi identities lead to a conservation equation for the total energy density. The two 
substances are separately conserved only if we assume that the two fluids does not interact or interact 
very weakly. This is a quite reasonable assumption since popular dark matter candidates (such as 
neutrinos and WIMPs) are indeed very weakly interacting particles. 
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As we know, observations tell us that we are living in a phase of accelerated expansion 
so that today the leading pressure must be negative. Let us first consider the case 7 > 
so that the pressure is positive (almost vanishing) in the limit pvdw << Pc and negative 
in the opposite limit pvdw >> Pc- This simple consideration suggests us that models 
with 7 > could have some chance to reproduce the observed data provided that the 
present day energy density of the dark matter is much larger than the Van der Waals 
critical energy density, i.e. it should be r] » 1 where hereon all quantities with a 
subscript are evaluated today (i.e. at z = 0). On the other hand, models with 7 < 
have a negative pressure in the limit 77 << 1 and a positive (and large) one in the limit 
77 >> 1. We may thus conclude that, in order to give an accelerated expansion today, 
t)q must be much smaller than 1 for models with 7 < 0. It is worth noting that these 
conclusions have been obtained without the need to solve the Friedmann equations. 

Data on the CMBR anisotropy spectrum strongly suggest that the universe is 
spatially flat today [2J EJ El so that we will assume k = in the rest of the paper. In 
this case, evaluating Eq.([TJ) at z — 0, we get : 

il c r] + £l b0 = 1 — > fi c = - — (8) 



with : 



^c = Pc/ Pcrit (9) 



with f^ )0 = Pb(z = 0)/ p C rit and p cr i t = 3Hq/8tiG the critical density of the universe. 
It is worthwhile to note that Eq.flEJ) is a consequence of the assumption that there is 
only a single fluid other than baryons filling the spatially flat universe and playing the 
role of both dark matter and dark energy. From this point of view, Van der Waals 
quintessence may be considered in the framework of unified dark energy models such 
as the Chaplygin gas ^3] and the Hobbit models jTHj. There is, however, a significant 
difference. The Van der Waals equation of state does not interpolate between a dust like 
pressure and a constant negative one as pvaw evolves. Indeed, for models with 7 > 
the pressure becomes quite small for 77 << 1, but it is not zero. On the other hand, 
when 7 < 0, the pressure becomes quite small and negative, but still remains different 
from zero. 

It is worth noting that it is possible to start narrowing the class of models to 
explore without the need to explicitly solve the Friedmann equations. To this aim, let 
us consider the deceleration parameter : 

ad a 1 13 p VdW 
a 1 a H l 2 2 p b + p VdW 

where we have combined Eqs.(JTJ) and (J2J. Inserting Eq.fl7J) into Eq. ffTUj) . evaluating the 
result at the present day and solving with respect to 7, we get : 

8(l-2g )(77o-3) 

7 9(l-a, )(3r/ 2 -9r ? o + 8) ' 1 ' 
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Since f2& >0 < 1, the denominator is always positive so that 7 is consistently defined by 
Eq. ffTTj) . We may thus use q as a parameter instead of 7 and characterize the models 
by the three parameters (qo, log 7/0, fib,o) where we use the logarithm of r/ instead of r] 
itself for reasons that will be clear later. It is worth noting that using q instead of 
7 allows to select immediately nowadays accelerating models. Moreover, it is easier to 
choose a range for q than for 7 that is a priori completely unknown. Actually, there 
is a third parameter that has to be given to completely assign the model and study its 
dynamics. This is the Hubble constant Hq which, however, we will consider as a known 
quantity. We will explain later why this assumption is necessary and why, nonetheless, 
it does not introduce any loss of generality. 

Finally, let us observe that there is another interesting quantity that may be 
evaluated without the need of solving the Friedmann equations. By definition, the 
barotropic factor w(z) of the Van der Waals fluid is given as : 

/ \ Pvdw 37 9 

w(z) = = — - -777O) . 

Pvdw 3 - r](z) 8 

Using Eq. (fTTj) and evaluating the result at the present day, we get : 

= "3(13?^) ' < 12 > 

For accelerating models (go < 0), wo takes on negative values and may also be lower 
than -1 thus suggesting that phantom models (that are indeed characterized by negative 
pressure fluids with wo < —1) could be somewhat confused with Van der Waals 
quintessence. 



3. The dynamics of the universe 



Having described the general features of the model, let us now determine the dynamics 
of the universe, i.e. let us investigate how the relevant physical quantities (scale factor, 
energy densities and Hubble parameter) evolve. As a preliminary step, let us remember 
that, solving the continuity equations for baryons, we get : 

n b ( z ) = n bfi (i + z f (13) 

with z — l/a — 1 the redshift. To solve the continuity equation for the Van der Waals 
dark matter, it is convenient to change variable from the cosmic time t to the redshift 
z and to use the dimensionless energy density 77(2). We thus get : 

dq = 377 [97r/ 2 -(277 + 8)r/ + 24(l + 7)] 

dz 8(l + z)(3-77) ' 1 ' 

This is a first order nonlinear differential equation for rj which may be numerically solved 
provided that the three parameters (q , log 770, S\ ) are given. It is worth stressing that, 

% We drop the subscript VdW because there is no possibility of confusion with the barotropic factor of 
the baryons that is identically zero. 
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Figure 1. The evolution against the redshift of the Van der Waals density parameter 
setting (gchlog^o) = (—0.38,-0.11) and three different values of Slbfih 2 , namely 
Slbflh 2 = (short dashed), n bfl h 2 = 0.0214 (solid) and n bfi h 2 = 0.0428 (long dashed). 
We use h = 0.664 to get f2f, i0 from Q bfi h 2 . 

although not explicitly present, the parameter fl^fl enters Eq. (jl4|) through 7. Moreover, 
being Eq. (f"H|) nonlinear, a sort of butterfly effect takes place with the result that also 
small changes in the baryons content lead to significantly different evolutions of the the 
energy density of the Van der Waals dark matter with the redshift z. This is clearly 
shown in Fig. [T] where we plot Qvdw( z ) = ^cV( z ) f° r a particular choice of (go, log 770) 
and three different values of VL^^h 2 . Comparing the short dashed line (corresponding to 
a model with no baryons) with the other two, we may safely conclude that, although the 
present day baryon density parameter is quite small (ttb,o ~ 0.04), it can not be neglected 
without introducing a severe bias in the past (z > 1) evolution of the Van der Waals 
energy density + . This is also evident considering the the Hubble parameter. While, for 
z < 1, there is essentially no difference at all among models having the same values 
of (go,logr?o) but different baryon content, the dependence on flb,o is more and more 
important going back in time, i.e. to higher redshifts. In particular, we note that, for a 
given z, the higher is flb,o, the larger is H(z), i.e. models with higher baryons content 
evolve faster. It is worth noting that these results do not depend on the values of the 

+ Note that we are defining Qvdw( z ) a s the ratio between the Van der Waals energy density and the 
present day critical density p cr i t — 3H 2 /8irG. It is thus not surprising that Qvdw(z) gets larger than 
1 even if the universe is assumed to be spatially flat. 
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Figure 2. The evolution against the redshift of the barotropic factor w for the same 
models considered in Fig.^ 



parameters (go, log ^o) that only determine when the effect of Qb.o becomes important. 

A striking feature of Fig.[T]is the flattening of flvdw(z) towards a constant value 
suggesting that the Van der Waals fluid mimics a cosmological constant term in the past. 
It is worth stressing that this behaviour does not depend on the particular 
choice of the parameters (qo, logrjo), but is rahter a characteristic feature 
of the Van der Waals fluid. However, the lower is the value of log 770, the 
higher is the redshift z when Qvdw{ z ) becomes constant. Although quite 
peculiar, this result may be easily explained considering the evolution of the 
barotropic factor w(z). As an example, Fig.[2]shows w(z) for the same models 
considered in Fig.[Tl First, we note that the effect of Q b0 is negligible (for 
realistic values of this quantity) in agreement with what is predicted for wq 
by Eq. (|12|) . Moreover, w(z) quickly converges to % = -1 so that the Van 
der Waals fluid behaves as the cosmological constant in the past. Changing 
the values of the parameters (go,log^o) have the only effect of increasing or 
decreasing the rate of convergence, but not the asymptotic limit of w(z). In 
particular, implicitly defining z\ as w(z\) = — 1, it is possible to show that z\ 
is a decreasing function of log 770, while it depends only weakly on q . Since 
Qvdw(z) becomes constant for z > z\, this explains why the energy density 
of the Van der Waals fluid becomes flat in Fig.[T] for higher values of z as 
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Figure 3. The evolution of the scale factor a with the scaled cosmic time t norrn = t/to 
for the same models considered in Fig.^ 

log 7/0 decreases. We may thus conclude that the Van der Waals fluid behaves 
as the cosmological constant at high redshift whatever are the values of the 
parameters (go, logr/o, fi&,o). This important result will have profound impact on the 
evaluation of the positions of the peaks of the CMBR anisotropy spectrum as we will 
discuss later in Sect. 6. 

Having determined H(z), it is straightforward to estimate the age of the universe 
at redshift z as : 

t(z) = / r — — dC = t H - —^-dC (15) 

1 } J z (i + c)#(C) h (i + C)^(C) 

where tu = c/H = 9.78/2T 1 Gyr is the Hubble time (being h the Hubble constant in 
units of 100 km s _1 Mpc" 1 ) and E{z) = H(z)/H . The present age of the universe is 
then to = t(z = 0). The evolution of the scale factor a as function of cosmic time t may 
be obtained by numerically inverting Eq. ()15|) and remembering that a = (1 + z)" 1 . The 
result is shown in Fig.Elfor the same models considered above. It is worth noting that 
a(t) does not enter anyone of the tests we will perform later that are only dependent 
on H(z). Nonetheless, it is interesting to look at the scale factor to get a feeling of how 
the universe evolves with time. 

It turns out that neglecting baryons leads to a strong error in the determination of 
the scale factor and thus on the age of the universe. The butterfly effect here works quite 
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hard. Actually, it is worth noting that the shape of a(t) does not depend strongly on 
the exact value of o as could be inferred comparing the solid and long dashed lines in 
Fig. El Doubling f2& iQ does not change significantly a(t/t ) so that we may conclude that 
what is important is to take into account the presence of baryons, but not the precise 
value of their density parameter (unless we take completely unrealistic values). 

4. The dimensionless coordinate distance 

A whatever model that aims at describing the evolution of the universe must be able 
to reproduce what is indeed observed. It is thus mandatory to test the viability 
of the proposed Van der Waals quintessence by contrasting and comparing it to the 
astrophysical data available up to now. This is also a powerful tool to constrain the 
model parameters thus paving the way towards a complete characterization of the model. 
As a first step, we may resort to the Hubble diagram of SNela, that is the plot of the 
distance modulus as function of the redshift z. However, we prefer here to follow a 
very similar, but more general approach considering as cosmo logical observable the 
dimensionless coordinate distance defined as : 



with E(z) defined above. For completeness, we remember that y is related to the usual 
luminosity distance (which is the quantity measured through the SNela distance 
modulus) as follows : 



It is worth noting that y(z) does not depend explicitly on H so that it is now clear why 
we are confident that our above choice for Hq does not alter the main result. Actually, 
H enters in the estimate of y bs(zi), the observed dimensionless coordinate distance to 
an object at redshift z t . Daly & Djorgovski ^7j have determined y(z) for the SNela in 
the Gold dataset of Riess et al. |3| which represents the most updated and homogeneous 
SNela sample today available. Since SNela allows to estimate rather than y, a value 
of H has to be set. Fitting the Hubble law to a large set of low redshift [z < 0.1) 
SNela, Daly & Djorgovski J7j have determined : 



We thus set h = 0.664 in order to be consistent with Daly & Djorgovski jTTj, but we have 
checked that varying h in the 68% CL quoted above does not alter the main results*. 

To increase the sample, Daly & Djorgovski added 20 further points on the y(z) 
diagram using a technique based on the angular dimension of radiogalaxies [T71 ITS] . 
Both this sample and the 157 SNela contained in the Riess et al. compilation span the 

* Note that the value we are using is consistent also with Hq = 72 ± 8 km s _1 Mpc -1 given by the 
HST Key project based on the local distance ladder. 




(16) 



D L = — (l + z)y{z) . 



o 



H = 66.4 ± 0.8 km s" 1 Mpc 



-i 
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redshift range (0.1, 1.8) so that it is possible to detect eventual systematic deviations of 
one tracer from another. None of such trends have been detected so that the full sample 
may be used without introducing spurious problematic features in the y(z) diagram. 
To determine the best fit parameters, we define the following merit function : 



where the observed quantities [z^y^ai) are given in Tables 1 and 2 of Ref. [T7j. 
Note that, although similar to the usal reduced x 2 introduced in statistics, the \ 2 
defined above is not forced to be 1 for the best fit model since the uncertainties Oi 
are not Gaussian distributed, but take care of both statistical errors and systematic 
uncertainties. Nonetheless, it is possible to compare different couples of model 
parameters on the base of the x 2 value. To determine constraints on the model 
parameters, we first define the marginalized likelihood functions : 



After having normalized to 1 at maximum, the 68% CL (95% CL) on a parameter pi is 
obtained by solving for Ci(pi) = exp (—0.5) (A(Pi) = ex P ( — 2)). 

Having set the Hubble constant as above, we are left with a three dimensional 
space to explore defined by the model parameters (go, logr/o, fi&,o)- Actually, we have 
seen that the main dynamical quantities we are interested in depend only weakly on 
Qb,o in the redshift range that is proven by the available astrophysical data. Moreover, 
flbfi is severely constrained by theoretical models of nucleosynthesis and by the observed 
abundance of light elements. Based on these considerations, Kirkman et al. |20j have 
estimated : 



Neglecting the small error, we thus set flb,oh 2 = 0.0214 and use the above quoted value 
of h to get our estimate of f^o- We have therefore only two parameters to constrain, 
namely go an d log?7o. Note that we use log 770 instead of t] itself since the former is 
easier to handle in numerical codes. As regard the range for go, since we are interested 
in accelerating models only, we set g = as upper limit. A lower limit may be obtained 
by the following argument. Let us insert Eq. lfTTf into Eq. flHJ) and evaluate the result 
at z = 0. We get : 




(17) 




(18) 




n bfi h 2 = 0.0214 ± 0.0020 . 




2=0 



2(1 + g ) - 3Q b>l 
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Figure 4. Observed data and theoretical curve of the dimensionless coordinate 
distance for the best fit Van der Waals model defined by (go, log ^o) = (—0.38, —0.11) 
with (h,Q bfi ) = (0.664,0.0214). 

It is reasonable to impose that the energy density is a decreasing function of cosmic 
time so that it is always dr]/dz > 0. In order to fulfill this condition at z = 0, we must 
impose g > (3£\ — 2)/2 — —1. Eq.(fT2|) makes it possible to see that the above lower 
limit on go excludes all models with wq < —1. While theoretically motivated by the 
need of not violating the energy conditions, such a choice seems to be disfavoured by 
the observations. Actually, fitting models with constant equation of state to the SNela 
Hubble diagram and the CMBR anisotropy spectrum gives wo = — 1. 43^38 121] • It is 
worth noting, however, that this estimate is strongly model dependent so that it does by 
no means imply that wq should be in that range also for the Van der Waals quintessence 
scenario we are investigating. Therefore, we have decided to use the above quoted lower 
limit for q although this excludes models with wo < —1. Nonetheless, we have checked 
that the constraints on go from observations are unaltered if we decrease the lower limit 
up to g = —2 thus entering the realm of w < — 1 models. 

There are no physical motivations to select a plausible range for log^o so that we 
have resorted to a trial and error procedure trying to not exclude any interesting regions. 
We thus end up with the following range for the model parameters : 

(3Q 6 , - 2)/2 < g < , -3.0 < logr/ < 0.0 . (19) 

By running (g ,log^o) m the range defined above, we find the best fit parameters: 
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(g , log vo) = (-0.38,-0.21) 

giving x 2 — 1-26. The fit is quite successful as can be seen in Fig.Hl while Fig. El shows 
the normalized likelihood functions from which we derive : 

- 0.52 < g < -0.33 , -0.42 < logr/ < -0.11 at 68% CL , 

(20) 

-0.58 < g < -0.28 , -1.38 < log rj < -0.05 at 95% CL . 

Note that the 95% upper and lower limits on go are significantly far away from those of 
the range (fT9")) which is a reassuring evidence suggesting that we have not excluded any 
observationally interesting region of the parameter space. While this is also true for the 
95% lower limit on log 770, the corresponding upper limit is quite near to the upper edge 
in Eq. ([T9|) which could be somewhat worrisome. However, the shape of the likelihood 
function for log r/ suggests that this is not a problem since the excluded region is outside 
the 99% CL. As a further consistency check, we have also repeated the fit adopting an 
higher upper limit on log 770, but the constraints (20) turn out to be unaltered. 

The first remarkable result of the fit is the exclusion of models with positive values 
of logr/o, i.e. models having a nowaday energy density of the Van der Waals fluid larger 
than the Van der Waals critical density (r/ > 1) are not compatible with the observed 
data. Moreover, the test is not very efficient to constrain the model parameters since 
only a minor part of the range for go is cut away, while logr/o is allowed to vary over 
almost an order of magnitude. Actually, this is not surprising. In the redshift range 
probed by our data (z < 1.8), the dimensionless coordinate distance for two models with 
very different values of log rj approximately coincide within the observational errors so 
that, to discriminate among distinct models, we need either more precise data in the 
redshift range (1.0, 2.0) or extending the dataset to z up to ~ 3. While it is likely that 
forthcoming satellite experiments such as the planned SNAP mission [22] will furnish 
more precise measurements of y in the redshift range (1.0, 2.0), going to higher redshift 
will need a tracer other than SNela. Good candidates in this sense are compact radio 
sources {23J since they can be detected up to z ~ 4, but there are still some problems 
related to the evolution with redshift of their physical properties. 

5. The age of the universe 

The present day age of the universe to is a powerful tool to constrain the parameters of 
a given cosmological model. Let us thus investigate what can be learned by applying 
this test to the Van der Waals scenario we are considering. 

First, we need an estimate of to from observational data. In Ref. [Zj, Rebolo et al. 
have performed a detailed combined analysis of the WMAP and VSA data on the CMBR 
anisotropy spectrum and SDSS galaxy clustering thus obtaining t = 14. 4^^ Gyr 
at 68% confidence limit. A more precise determination has been obtained by Seljak 
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Figure 5. Marginalized likelihood functions (normalized at 1 at maximum) for the 
two model parameters qo (left panel) and log ryo (right panel) . 



et al. [21] who have fitted the ACDM model to a combined dataset comprising the 
CMBR anisotropy spectrum, the galaxy power spectrum, the SNela Hubble diagram, 
the dependence of the galaxy bias on mass and the Lya clouds power spectrum. As 
a result, they get t = 13.6 ± 0.19 Gyr at 68% CL. Actually, both these estimates are 
model dependent since the authors assume a priori a background cosmological model 
and then determine to from the best fit parameters of that model. In order to avoid any 
systematic bias, we prefer to use a model independent estimate although this leads to 
enlarge the error bars. To this end, we resort to age estimates of globular clusters and, 
following Krauss > we retain as viable all models such that : 

to e (10.2, 16.0) Gyr . (21) 

Note that the range quoted above is consistent both with the model dependent estimates 
of Rebolo et al. and Seljak et al. and with t > 12.5±3.5 Gyr from radioisotopes studies 
[2*o] so that we are confident not to be excluding interesting models. 

Given (g ,log?7o), the predicted age of the universe for the corresponding Van der 
Waals model may be estimated setting z = in Eq.(fT5j). Fig. shows the age contours 
in the (go, log ^o) plane for the model with (h,Qb,o) set as discussed above. Note that 
we have greatly enlarged the range for log 7/0 in order to explore with much detail this 
poorly constrained parameter. The results of the age test are somewhat surprising. 
First, positive values of log i] are strongly excluded since they give rise to values of to 
of the order of hundreds of Gyr in striking disagreement with Eq.(|21|). Furthermore, 
the age test also cuts away a large part of the range for q accepting only models with 
qo > —0.29. Roughly, we can summarize the results of this test giving the following 
constraints : 

- 0.29 < g < , -6.0 < logr/ < -0.5 . (22) 

Note that this is actually a little bit larger than the region delimited by the dashed line 
in Fig. H since we have approximated it as a rectangular one, while it is not. However, 
this simplification will not alter our main results. Note that, for the model parameters 
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Figure 6. Contour plots of equal to in the plane (goilog??o)- The solid lines refer to 
values from 14.5 to 16.5 Gyr (in steps of 1 Gyr) from right to left. The dashed line is 
the observed upper limit to — 16.0 Gyr. 

in the range (J22J) . pvdW,o — lPvdw,o with 7 < so that the Van der Waals fluid behaves 
essentially as a perfect fluid with negative pressure. 

Comparing Eq. fl22R with the constraints (20) from the fit to the 
dimensionless coordinate distance shows that the two ranges have a good 
overlap in log 770 a t 95% CL, but the constraints on q are marginally consistent 
only at the 95% CL. Moreover, all the allowed models predict values of to close the 
upper end in Eq. (j21j) . One should thus conclude that, in order to reconcile Van der 
Waals quintessence with both the SNela and radiogalaxies distance data and the age of 
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the universe, a certain degree of fine tuning of the parameters is needed. Although this 
is not a problem in its own, it is worth noting that the Van der Waals equation of state 
is only an approximation to the actual unknown equation of state of the dark matter 
component which is likely to be valid only over a limited range of the evolution history of 
the universe. Unfortunately, a detailed knowledge of the dark matter thermodynamical 
properties and of their possible dependence on z is completely lacking so that it is 
impossible to say what is the redshift range over which Eq.((7j) holds. This systematic 
error does not affect the fit to the dimensionless coordinate distance since this test 
only probes a limited and nearby redshift range, but could have a profound impact 
on the estimate of to since this latter depends on the full evolutionary history. As a 
consequence, we do not consider as a serious shortcoming the only marginal agreement 
of predicted and estimated age of the universe for the Van der Waals models selected 
by the fit to the SNela and radiogalaxies dimensionless coordinate distances data. 

6. The peaks of the CMBR and some remarks on structure formation 

The cosmological model we have discussed insofar is made out of matter only, but we 
have made a clear separation between the baryons and dark matter. While the former 
is still described as dust matter, the equation of state of the latter is the Van der Waals 
one so that the properties of the dark matter in this scenario are radically different from 
that of the standard cold dark matter (CDM). Indeed, while for CDM it is w = at all 
z, for the Van der Waals dark matter w depends on z and quickly approaches the value 
w = — 1 as shown in Fig.El It thus makes sense to ask how structure formation evolves 
in the scheme we are proposing. This is a quite complicated task and will be addressed 
in detail in a forthcoming paper. Nonetheless, here we give some qualitative comments 
to illustrate the subtleties of this topic. 

As a first remark, let us remember that structure formation may efficiently take 
place only during a decelerating phase of the universe evolution. Moreover, the SNela 
Hubble diagram shows some evidences of a transition from acceleration to deceleration 
although the estimate of the transition redshift zt (defined so that q(zr) = 0) are 
quite model dependent. For instance, Riess et al. 5 J obtain z T = 0.46 ± 0.13 by 
using the ansatz q(z) = q + (dq/dz) z=0 z, while the detailed analysis of the ACDM 
model performed by Seljak et al. [24J gives 0.52 < zt < 0.91. In Fig.0 we report the 
deceleration parameter for the best fit model parameters (go,log^o) an d three different 
baryon contents. For this particular model, q(z) takes positive values only in the redshift 
range (~ 0.4, ~ 0.7) so that the universe is accelerating also during the period when 
structure formation takes place in the concordance ACDM model. The situation is 
partially ameliorated considering models with larger values of log?7o since the redshift 
range over which q(z) is positive enlarges and shifts to higher z. However, one should also 
take into account the effect of Qbfih 2 since this could significantly impact the behaviour 
of the deceleration parameter. 

Likely, the most worrisome feature of the Van der Waals quintessence scenario is 
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Figure 7. Deceleration parameter q as function of the redshift z for the best fit 
model with (qoAogVo) — (—0.38, —0.21) and three different values of fl^oh 2 , namely 
Slbflh 2 = (short dashed), n bfl h 2 = 0.0214 (solid) and Q. bfl h 2 = 0.0428 (long dashed). 
We set h = 0.664 to convert from Qb,oh 2 to Slb,o- 



the barotropic factor w going to -1 in the past since this makes the Van der Waals 
fluid never behaving as standard cold dark matter which could suggest a completely 
unrealistic growth of perturbations in our scenario. However, such a conclusion is at 
least premature. Let us consider the ACDM model that correctly describes the evolution 
of structures as we observe it. There are two fluids, standard cold dark matter with p = 
and the cosmological constant with p = —p, and during the structure formation epoch 
the CDM energy density dominates over that of the cosmological constant. Actually, a 
very similar situation takes place in our model. Indeed, we have checked that, in the 
redshift range where presumably structure formation takes place, the baryons energy 
density dominates over that of the Van der Waals fluid that, in this period, is very 
well approximated by a cosmological constant - like term for all values of the model 
parameters. Therefore, in the far past our model is formally equivalent to the ACDM 
model with the baryons and the Van der Waals dark matter playing the roles of CDM 
and A respectively. This nice result suggests that structure formation could evolve in a 
very similar way, but a detailed investigation is needed to draw a definitive answer. 

As a first preliminary investigation, we may compute the positions of the first three 
peaks in the CMBR anisotropy spectrum using the procedure detailed in Refs. [27| 12%]. 
According to this prescription, in a flat universe made out of a matter term and a scalar 
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field -like fluid, the position of the m-th peak is given by: 



lm = Ia{™> - <P - 



(23) 



with I a the acoustic scale, (p the overall peak shift and 5ip m the relative shift of the 
m-th peak with respect to the first. While <p and 5<p m are given by the approximated 
formulae in Ref. j2H], the acoustic scale for flat universes may be evaluated as [27| : 
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In Eqs. (J24)l and (j2SJ), «Zs = (1 + 1 with z/ s the redshift of last scattering that may 
be computed as [2H] : 
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with ((71,^2) given in Ref. 
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and (j2H) are 

(27) 
(28) 
(29) 



where r — J a~ x dt is the conformal time, and p r are the energy densities of the 
baryons and radiation respectively, w(z) and = p^ / p cr it(z) are the barotropic factor 
and the density parameter)) of the dark energy fluid. In order to use Eqs. ()24j) - (|29j) . we 
first remind that the only matter term in our model is the baryonic term so that we 
have to replace everywhere Qm with 0&. On the other hand, the role of the scalar field 
fluid is played by the Van der Waals fluid so that all the quantities with the subscript </> 
have now to be evaluated using the Van der Waals energy density. Finally, the Hubble 
parameter is given by Eq.(J2J) with pvdw( z ) obtained by solving Eq.([14|). while the present 
day value of the radiation density parameter is set as Qq = 9.89 x 10~ 5 . Note also that 
one has to choose a value for the index n of the spectrum of primordial fluctuations 
which we set as n — 1. 

jj Note that here we are using p cr u{z) — 3H 2 (z)/8irG rather than its present day value which has been 
adopted to get Fig.^ As a consequence, in this section, we accordingly redefine Qvdw{z) as the ratio 
between pv<w{z) and p cr u(z). 
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Figure 8. Constraints in the (qo,\ogi]o) plane for models with (h,Qb,oh 2 ) — 
(0.664,0.0214). The region to the right of the solid line individuates models that 
are consistent with l^ oom within 3<r. Models with (go,log?7o) belonging to the region 
of the parameter space to the right of the short and long dashed lines agree within 
5cr with I^map anc j jWMAP reS p ec tively. We only report the region of the parameter 
space delimited by the 95% CL constraints in Eq.(20). 



The position of the first two peaks in the CMBR anisotropy spectrum ha been 
determined with great accuracy by the WMAP measurements giving [U] : 

jWmap _ 220.1 i 0.08 , /f MAP = 546±10 , (30) 
while the position of the third peak is more uncertain and may be estimated as [2] : 

Fig. |H1 shows the region of the parameter space (q , log rj ) that are consistent with the 
bounds from the position of the peaks. Note that the model predicted l m are affected by 
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unknown errors due to both the approximated nature of the procedure adoptedff and to 
some subtleties we discuss later. To qualitatively take into account these systematics, we 
generously enlarge the error bars on the measured l m considering the 5a (3a) confidence 
ranges for WMAP (BOOMERanG) observed quantities. Considering the most stringent 
cut (that on li), Fig. |H1 shows that it is indeed possible to find out models that are in 
agreement with both the fit to the dimensionless coordinate distances and the position 
of the first three peaks. However, some caution is needed. 

The method adopted to compute l m implicitly assume that the standard theory 
of perturbation may be applied to the Van der Waals quintessence scenario we are 
discussing. In this case, the sound speed c s entering the growth of perturbations is 
usually defined as c 2 s /c 2 — p/p with c the speed of light. Using Eqs.(HJ), we easily get : 

■c.\'_ 9 7 (4 -,)(1 -i,) 2 



V c /( 



cJ 4(3 - ?7) 2 

that at z = reduces to : 

2 6(77o-4)(77q-1) 2 k;o 
ch (r ]o -3)(3 V 2 -9r ]o + 8) 1 ' 

having used Eq.fjllj) . This quantity is not well behaved since, for the model parameters 
in the range determined above, may become negative giving rise to a formally imaginary 
sound speed. Moreover, during the evolution of the universe, c 2 s can take negative values 
over a large redshift range. From this point of view, Van der Waals quintessence is similar 
to other unified dark energy models. In particular, Sandvik et al. [3D] have considered 
the particular case of the generalized Chaplygin gas and shown that it gives rise to 
oscillations or exponential blowup in the dark matter power spectrum inconsistent with 
observations. Although only the generalized Chaplygin gas has been investigated, they 
argue that similar problems also take place for every unified dark energy model because 
of these models having a negative c 2 s . This could be an evidence against the Van der 
Waals scenario we are proposing. Actually, things get different. First, one has to remind 
that p/ p is the adiabatic sound speed. If the fluid fluctuations were adiabatic, then the 
pressure perturbation Sp corresponding to a perturbation Sp in the energy density should 
be |2H : 

Sp = (p/p)Sp (34) 

so that if p/p becomes negative or singular the perturbations go unstable. This is the 
case for our model, but this problem may be avoided by resorting to non adiabatic 
pressure perturbation. In this case, the relation between Sp and Sp is generalized as: 

Sp = (p/ p)Sp + pT (35) 

ft To test the validity of the method adopted to compute l m , one should compare the predicted values 
with what is obtained running a code that evaluates the full spectrum of perturbations. This test has 
been indeed performed with success only for scalar field dark energy with different choices of the self 
interaction potential. Although it is likely that the method works with the same goodness also in our 
case given the similarity of our model with the cosmological constant in the prerecombination epoch, 
a detailed comparison is still needed. 
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with T the entropy contribution related to the anisotropic scalar stress ir L . Using the 
freedom in the choice of tt^, it is then possible to define an effective sound speed that 
determines the scale under which the fluid energy is effectively smooth through the 
sound horizon. Well above this scale, stress gradients are negligible and the standard 
theory of perturbations applies. It is thus interesting to investigate how it is possible 
to implement such a mechanism for the Van der Waals quintessence and then study the 
problem of structure formation in this generalized scenario. We will address this topic 
in a forthcoming paper [32J. 

7. Conclusions 

The increasing bulk of astrophysical data accumulated in recent years has delineated a 
new standard cosmological paradigm. According to this picture, the universe is spatially 
flat and driven by an unknown form of dark energy leading to an accelerated expansion. 
Soon after the establishment of such a scenario, the hunt for candidates to the dark 
energy throne has started leading to the proposal of a plethora of mechanisms ranging 
from the old cosmological constant to various scalar field quintessence, to modification 
of Friedmann equations (motivated by extradimensions and braneworld theories) and 
higher order geometrical terms in the gravity Lagrangian [UHl EH EH] • Although being 
completely different in their dynamical properties and underlying physics, they all share 
the ability of well fitting the same set of data. 

Furthermore, it is worth noting that the dark side of the universe is not only 
populated by dark energy, but also by dark matter whose nature is still far to be 
understood. It is thus tempting to ask whether these two ingredients are indeed different 
substances or two aspects of a single fluid whose properties are different from what we 
usually expect. To this regard, it makes sense wondering if the dust approximation 
used for matter is the right one to describe the equation of state of the dark matter. 
In a tentative to explain the observable quantities of our universe with the minimal 
number of ingredients, here we have investigated the possibility that dark matter and 
dark energy are actually a single fluid whose equation of state is that of a Van der Waals 
gas. In Van der Waals quintessence scenario, there is a single fluid whose equation of 
state comes directly from classical thermodynamics since the perfect gas approximation 
cannot be used to describe phase transitions which occur during the evolution of the 
universe. Although their contribution to the energy budget is nowaday subdominant, 
we also include baryons in the model since, because of the nonlinear character of the 
dynamical equations, they play an important role to determine how the main quantities 
evolve with the redshift. 

Any theory, as elegant and motivated it can be, is meaningless if it is unable to 
give a coherent description of the universe as it is observed. That is why we have tested 
Van der Waals quintessence against the dimensionless coordinate distance to SNela and 
radio galaxies. This allows to narrow the parameter space of the model considering the 
present day values of the deceleration parameter and of the ratio between the critical 
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density of the universe and the Van der Waals critical density. We find : 

- 0.52 < g < -0.33 , -0.42 < logr/ < -0.11 at 68% CL , 

-0.58 < g < -0.28 , -1.38 < log r] < -0.05 at 95% CL , 

with (g ,log?7o) = (—0.38,-0.21) as best fit parameters. The constraints on logr/o are 
weak because of a serious degeneracy among the model parameters. Indeed, different 
models of this class predicts values of y(z) which are in agreement with each other 
within the observational errors. However, lowering the uncertainties on y(z) in the 
redshift range (1.0, 2.0) or extending observations to higher redshifts (up to z ~ 3 — 4) 
will allow to break this degeneracy. Both these possibilities are likely to become realistic 
in a near future thanks to the next -to -come satellite and ground based experiments. 

Motivated by this encouraging result, we have compared the predicted age of 
the universe with that estimated from globular clusters in a model independent way. 
Unfortunately, it turns out that only a quite small region of the parameter space 
delimited by the above constraints is allowed by this test. Moreover, t turns out to 
be systematically high so that one could argue against the proposed model. Actually, 
this conclusion is premature since our ignorance of the thermo dynamical properties of 
the dark fluid makes it impossible to determine the redshift range over which Eq.((7|) 
holds. Indeed, during the evolution of the universe, the dark fluid may undergo phase 
transitions thus changing its equation of state in a unpredictable way. Since to depends 
on the full evolutionary history of the universe, using Eq.(j?j) at all z may induce a 
systematic error if the relation p = p(p) is altered. Given these theoretical uncertainties, 
overemphasizing the result of the age test should be avoided. 

As a further test, we have also evaluated the position of the first three peaks 
of the CMBR anisotropy spectrum using the analytical approximation developed in 
Refs. [27J I2H]- Comparing with the measured l m , we have individuated a region of the 
parameter space (go, log 770) which is consistent with both the distance fit and the peaks 
position. Although very encouraging, this result should be confirmed by evaluating the 
full spectrum and not the peaks position only. To this aim, one has to first investigate 
in more detail how perturbations grow in the Van der Waals quintessence scenario. 
Because of the negative adiabatic sound speed, perturbations go unstable and one has 
to take into account also the role of anisotropic stresses and entropy production. 

We would like to conclude with a general comment. Van der Waals quintessence has 
turned out to be an interesting scenario for describing the late universe and seems able 
to solve the puzzle of dark energy without adding exotic fluids or arbitrary modifications 
of Friedmann equations. On the other hand, classical thermodynamics tells us that the 
Van der Waals equation of state is only an approximated description of a realistic fluid. 
In our opinion, before invoking the help of new physics, it is worth wondering whether 
classical physics could still suggest the way to shed light on the dark side of the universe. 
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